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Abstract 

We appeal to a complex q-Fourier transform as a generalization of the 
(real) one analyzed in [Milan J. Math. 76 (2008) 307]. By recourse to 
tempered ultra-distributions we are able to show that the q-Gaussian 
distribution can be nicely obtained, overcoming all troubles that afflict 
its real counterpart. 

KEYWORDS: q-Fourier transform, tempered ultradistributions, complex- 
plane generalization, one-to-one character. 
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1 Introduction 



Nonextensive statistical mechanics (NEXT) [TJ |2j [3J, a well known gener- 
alization of the Boltzmann-Gibbs (BG) one, is used in many scientific and 
technological endeavors. NEXT central concept is that of a nonadditive 
(though extensive [3]) entropic information measure characterized by the 
real index q (with q = 1 recovering the standard BG entropy). Applica- 
tions include cold atoms in dissipative optical lattices [5], dusty plasmas 
[BJ, trapped ions [7J, spin glasses [8j, turbulence in the heliosphere [9], self- 
organized criticality [10], high-energy experiments at LHC/CMS/CERN [TT] 
and RHIC/PHENIX/Brookhaven [12], low- dimensional dissipative maps [13], 
finance [14], galaxies [15], and Fokker-Planck equation's studies [16J, etc. 
NEXT can be advantageously expressed via q-generalizations of standard 
mathematical concepts [TT]. One can mention, for instance, the logarithm 
and exponential functions, addition and multiplication, Fourier transform 
(FT), the Central Limit |18J . plane waves, and the representation of the 
Dirac delta into plane waves [201 1211 1221 123] ■ 

Until recently, a generic analytical expression for the inverse q-FT for ar- 
bitrary functions and any value of q did not exist [23]. This situation 
was adequately remedied in [23], whose authors, by using tempered ultra- 
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distributions [121120], introduced a complex q- Fourier transform F(k, q) which 
exhibits nice properties and is one-to-one. In turn, this overcame a serious 
flaw of the original F q — definition, i.e., not being of the essential one-to-one 
nature Investigations of this kind and related questions are relevant 

for field theory and condensed matter physics, engineering (e.g., image and 
signal processing), and mathematical areas for which the standard FT and 
its inverse play important roles. 

In this work we focus attention on q-Gaussians, an essential tool of q-statistics 
[25] . that was not discussed in [25]. q-Gaussian behavior is often encountered 
in quite distinct settings [26]. In particular, one has to mention experimental 
scenarios in which data are gathered using a set-up that performs a nor- 
malization preprocessing. The ensuing normalized input, as recorded by the 
measurement device, will always be q-Gaussian distributed, if the incom- 
ing data exhibit elliptical symmetry, a rather common feature [26]. The 
q- Fourier transform of the q-Gaussian was discussed in [18j, but the corre- 
sponding treatment also presented the flaws above alluded to, a situation 
deserving further discussion that we tackle below. 
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2 Preliminaries 



It is necessary, before proceeding, to review materials developed in [25] (more 
details in the Appendix). So-called q-exponentials 

e q (x) = [1+(1-q)x]y (1 - q) , (2.1) 

are the hallmark of Tsallis's statistics [Tj, being generalizations of the ordi- 
nary exponential functions and coinciding with them for q = 1 . Here we will 
deal with complex q-exponentials, i.e., e q (ikx) for 1 < q < 2 with k a real 
number (see ([23]) 

e q (ikx) = [1 +i(1 - q)kx]T=q. (2.2) 

Our central tools are distributions, that is, linear functionals that map a set 
of conventional and well-behaved functions, called test functions, onto the 
set of real (complex) numbers. The best known example is perhaps Dirac's 
delta, the derivative of the Heaviside step function H 



H(x) = < 



1 for x > 

(2.3) 



for x < 0. 

Tempered distributions constitute a subset of the distributions-set for which 



the test functions are members of a special space called Schwartz' one S, 
a function-space in which its members possess derivatives that are rapidly 
decreasing. S exhibits a notable property: the Fourier transform is an au- 
tomorphism on S, a property that allows, by duality, to define the Fourier 
transform for elements in the dual space of S. This dual is the space of 
tempered distributions. In physics it is not uncommon to face functions that 
grow exponentially in space or time. In such circumstances Schwartz' space 
of tempered distributions is too restrictive. Instead, ultra-distributions sat- 
isfy that need [T9] , being continuous linear functionals defined on the space 
of entire functions rapidly decreasing on straight lines parallel to the real axis 

mi- 

An important fact about ultra-distributions is the following: a tempered dis- 
tribution is the cut of a tempered ultra- distribution. We are not speaking 
here of the "cut" of an analytic function (see Refs. [301 EI])- Accordingly, 
e q (ikx) is the cut along the real k tempered ultra-distribution [301 EI] > 

an essential fact for our present endeavor, 



E q (ikx) ={H(x)H[3(k)]-H(-x)H[-3(k)]}[l -q)kx]^f, (2.4) 

with H(x) the Heaviside's step function and 3(k) the imaginary part of the 
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complex number k. The relationship between e q (ikx) and E q (ikx) becomes 
more clear noting that, if f(k) is a tempered distribution and F(k) is the 
corresponding tempered ultra-distribution, then ([25 | [30 | I3T]) 

f(k) = F(k + iO) -F(k-iO). (2.5) 

Thus, (12. 5p leads to 

e q (ikx) = F_ q [i(k + iO)x] - E q [i(k - iO)x] (2.6) 
At this stage we introduce the set A[]^),oo, defined as 

An,2),oo = {f W/f (x) G A+ 2)>00 A f (x) G A^J, (2.7) 

where 

A n,2),oo = {f(x)/f(x){l +i(1 - q)kx[f(x)]t"-^}^ G C ] [R + ]A 

[f(x) >0;1 < q <2]}, (2.8) 

and 

A n,2),oo = {f(x)/f(x){1 - q)kx[f(x)]M}A G C'[R~]A 

[f(x) > 0;1 < q <2]}. (2.9) 
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With the help A and by recourse to (12. 4p together with the fact that for a 
given F(k, q) 



so that 



lim 

e->0+ 



dqS(q-1-e)F(k,q) = F(k) 



(2-10) 



f(x) = — cb dkF(k)e- ikx , 
2n 



(2.11) 



we can define a complex Umarov-Tsallis-Steinberg (UTS) q-Fourier trans- 
form (of f(x) G A[i ! 2) ; oo) in the following way 



F(k,q) = [H(q-l)-H(q-2)]x 



HCJ(k)] 



f(x){1 +i(1 -q)kx[f(x)] (q - 1] }^, dx- 



H[-3(k)] 



f(x){l +i(1 - q)kx[f(x)] { ^ ] }^ dx 



(2.12) 



Here q is a real variable such that 1 < q < 2. The cut along the real axis 
of this transform is the real UTS q-Fourier transform given in JTffl . f2^ (see 
I25tf for a simple application of this transform). Taking into account that for 
q = 1 the q-Fourier transform is the usual Fourier transform and using the 
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formula for the inversion of the complex Fourier transform straightforwardly 

leads to the inversion formula for 112.12]) . 

Consider 



i-2 



F(k) = lim 

e->0+ 



6(q-1 -e)F(k,q)dq, 



together with 



f (x) = T - 
2n 



dkF(k)e 



-ikx 



Since for q = 1 our equation 112.12]) is the complex Fourier transform 



F(k) = Hcr(k)] 



from $2.12]) we find 



r 



dxf(x)e lkx -H[-3(k)] 



dxf(x)e lkx , 



f(x) 



1 

2n 



o 



lim 
r L 1 



Akx dk. 



(2-13) 



F(k,q)5(q-1-e) dq 

Eqs. 112.12]) and A2.13]) solve the problem of inversion of the q-Fourier trans- 
form, which is of a one-to-one nature (see [27] for fixed q). Clearly, from 
f l23|) and (l2~lj) . on the real axis, one gets for fl2TT2|) and ( 12TT3|) 



F(k,q) = [H(q-1)-H(q-2)]x 
f(x){l +i(1 - q)kx[f{x)] (q - ]] }^ dx, 
for the real transform, and 



(2.14) 



f (x) = ^ 
2n 



lim 

e->0+ 



F(k,q)6(q-1-e) dq 



e~ ikx dk, 



(2.15) 



for its inverse. 



3 Series expansion of the q- Fourier transform 



Consider now the function 



e q [ikxf{x)^] = H(x,k,q) = {1 + (1 - q)ikx[f (x)]^ 1 }^, (3.1) 

that constitutes a generalization of e q [ikx], whose treatment was discussed 
above. Using the series expansions of the logarithm and the exponential 
function, we write 

ejikxflx)^ 1 ] ee {1 + (1 - qHkxtfM]^ 1 }^ = eT Vn{i+(i-q)iicx[fM]«-i} = 



e n=1 



1 £ (-^ n+1 (1 _q)n tikx) n e n(q-1)lnf(x) 

e n=1 



g H n=1 m=0 



oo co m— 1 

L L 2-^(ikx)"ln m [f(x)](q-1)"+— 1 
n=1 m=0 



{£ 

> ln=0 



(ikx) n - m+1 ln m [f(x)] 



(q-D n } 



(3.2) 
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Let g(x, k, n) be given by 



g(x,k,n) = ^ 1 _ 7 J (ikx) n - m+1 ln m [f(x)]. 



m=0 



m! 



Then, 



e q [ikxf(x) q_1 ] - q)ikx[f (x)]^ 1 } = H(x, k, q) = e 



L g(x 



or 



i -i Y. g(x,k,n) (q— 1 ) r 

eqtikxfM^ 1 ] = H(x,k,q) = e lkx e^ 



According to the exponential function's expansion we have 



... , Y_ qfx, k, n)(q — 1) n 

g n=1 v 



p=0 



so that 



e q [ikxf(x) q ] ] = h(x, k, q) = e 



ikx 



1 +^l(x,k,n) 



n=1 



where 



oo s— n+1 s— si — n+2 s— si— S2 s n _2~ 1 



Kx,k,n) = -XZ Z Z 

s=n si =1 S2=l s n _i=l 

g(x,k,si)g(x,k, s 2 ) • • • g(x,k, s n _i)x 



g(x,k,s-si -s 2 Sn^Hq-l) 5 
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Finally (see Section 2) we can write the q-Fourier transform of e q [ikxf (x) q ] ] 
in the fashion 



F(k,q) = [H(q-1)-H(q-2)]x 



H[3(k)] 



f(x)h(x,k,q) dx-H[-3(k)] 



f(x)h(x,k,q) dx } . (3.9) 



4 The q-Fourier transform of the q-Gaussian 

As stated in the Introduction, our purpose is to calculate the q-Fourier trans- 
form of the q-Gaussian. The calculation is too involved if one wishes to con- 
sider the pertinent expansions up to arbitrary order on q — 1 . We will content 
ourselves here with a first order approach. Accordingly 



h(x,k,q) = e ikx [1+g(x,k,1)(q-1)], 



(4.1) 



with 



g(x,k,1) = ^£+tkxln[f(x)], 



so that, up to first order, we have for the q-Fourier transform 



(4.2) 



Hp(k)] 



F(k,q) = [H(q-1)-H(q-2)]x 
f (ikx) 2 



1 + 



+ ikxln[f(x)] [> (q - 1) } f(x)e ikx dx 
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H[-3(k)] 



1 + + ikxln[f(x)] \ (q - 1) [> f(x)e ikx dx J> (4.3) 



i— 



Let G(k) and G(k, (3) be given by: 



G(k) = <^ H[a(k)] 



f(x)e ikx dx-H[-a(k)] 



f{x)e lkx dx } (4.4) 



G(k,|3) = jHp(k)] 
enabling us to write 



[f(x)]<V kx dx-H[-3(k)] 



[f(x)] p e ikx dx } , (4.5) 



F(k,q) = [H(q-1)-H(q-2)]x 



G(k) 



(q-1). (4.6) 



Let f (x) be the q-Gaussian 



f(x) = C,[1+(q'-1)ax 2 ]^7, 



(4.7) 



where: 



B 



1 _L_1 

2> q'-l 2 
Ci = 



Thus, using [28] we obtain 



G(k,q') =H[3(k)]C 



(4.8) 
(4.9) 



q 2 



[(q'-l)cx]^ 7 
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(1 — q')iak 



2-1 1 

i-q' 



X 



H 



2-q' 1 



ik 



N 



, q .-H(i-q'W f^-HO-q> 



ik 



H[-a(k)]c 



q 2 



H 



[(q'-1)a]i-< 
lk ' N 



T^-H(q'-U« 



(q'-l)ictk 
ik 



2-q 

1-q 



X 



7^ V(q'-D« 



(4.10) 



r f P+i-q 

7t ' ^ 1-q' 



G(k,q',(3) = H[a(k)]C p , 

2 [(q'-1)a]i-« 

L - L ^- \ V(1-q')a 



i-q' 



ik 



(1 — q')i<xk 
ik 



p+'-q 1 



i-q 



* — I 



1-q 



Vd-q> 



H[-3(k)]C p , 



H, 



r f P+i-q 

q' 7 — 

Z [(q'-l)a]^ 7 



2 



ik 



N, 



(q' — 1 )iak 
ik 



p+!-q 1 

1-q 



X 



X 



(4.11) 



where H and N are the Struve and Neumann functions, respectively. The 
q- Fourier transform of the q- Gaussian is now 



F(k,q,q') = [H(q-1)-H(q-2)]x 



G(k,q') + 



(q-1). (4.12) 



(5=1 



The cuts on the real axis of G(k, q') and G(k, q , (3) are 

2 



G(k,q ) = C 



71 



r f 2-q 



q 2 



[(q'-1)a]^7 
14 



(1 - q')ia(k + i0) 



2-q, i 
7^7 2 



X 



H 



i(k + iO) 



N 



C 



71 r (ff' 



q 2 

H 



[(q'-l)ct]^ 7 
i(k-iO) 
(q'-l)cx 



t(k + iO) 
&-H(l-q')« 

2 



+ 



2-q 1 
1-q' 2 



2-q 1 

2 



(q'-l)ia(k-iO) 

i(k-i9) 



x 



N 



G(k,q',(3) = C q p , , 



- r ( ^ 



[(q'-1)a]i-«< 
i(k + iO) 



-N, 



2 

(1 - q')ia(k + iO) 
t(k + tO) 



(4.13) 



1 

i-q' 



+ 



C p , 



p V i-q 



H 



q' ? -2-, 
i(k-iO) 



p+'-q I 



l-q 



f 2 



(q' - l)ia(k-iO) 

i(k-iO) 



x 



N, 



e±i-L_iv („> i \ ™ V(q'-i)a 



, (l H(q'-iK 

Now, the real q-Fourier transform takes the form 



x 



(4.14) 



G(k,q') + 



F(k,q,q') = [H(q-1)-H(q-2)]x 
k 2 3 2 /. , 3 3 / „. 



(q-1). (4-15) 



5 The q-Fourier transform of the q-Gaussian 



for fixed q 



In this section we will provide an alternative path to the computation of the 
q-Fourier transform, for a q-Gaussian, in the case of a fixed q-value. We start 
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with 



F(k,q) =Hp(k)] 



C q [l + (q-1)ax 2 ]^x 



i 

1 + (1 - q)ikx \ C q [1 + (q - 1 )ctx 2 ]^ >> (q - 1 |> 1_q dx 



H[-3(k)] 



C q [l +(q-1)ax 2 ] — x 



1 + (1 -q)ikx{C q [1 + (q-1)ax 2 ]— Wq - H " dx, (5.1) 



1 < q < 2. Simplifying terms we obtain 



F(k, q) = H[3(k)] 



C q [(q - 1 )ccx 2 + e-^(q - 1 )C^kx + 1]^ dx- 



H[-2f(k)] 



C q [(q -l)ax 2 + e^(q - 1)C q H kx + 1]^? dx. (5.2) 



Effecting the change of variables \/ iq — 1 )a x = y the q-Fourier transform 
adopts the appearance 

Hp(k)] 



F(k,q) 



^/[q^Tja. 
H[-a(k)] 



C f 



ky + 1 



dy- 



y 2 + e-q- 



q-i 



ky + 1 



i 

'-q 



dy. (5.3) 



Using now the result given in [33] (P!^ is the associated Legendre function) 



2^r(1 -2^t)(z 2 -1)i 

r(i-^)r(-n-v)r(v-^ + i) 
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(1 + 2tz + 1 2 ) H 1 -r 1 -^ dt, (5.4) 



and we can write: 



(l+2tz + t 2 f I t- 1 - v ^dt = r(-^)2-^ 1 (z 2 -1)2P^_ 1 (z), (5.5) 



where: 



q- 1 /q^T ii 



so that 



F(k, q ) = C q ^-V] 2-^ (y 2 k 2 + 1 ) * x 

v(q - i)a 

{H[3(k)]PVje-^yk) - H[-J(k)]e^P^^(e^yk)} , (5.7) 

which is the q- Fourier transform of the q- Gaussian on the complex plane for 
fixed q. 

Our next step is to evaluate the q-Fourier transform of the q-Gaussian on the 
real axis. This needs calculating the cut of (15. 7p along the real axis. Using 
we reach the equality 

(y 2 k 2 + 1)£ {Hp(k)]P^_Je^yk) - H[-3(k)]e l ^P^_Je^yk)} = 



H[-J(k)]e^ ( ff l] J (-ii, 1 + ; 1 - W ^— ^) =C(k,q), (5.8) 
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where F is the hypergeometric function. 

The cut c(k, q) along the real axis of C(k, q) is 



r , s (yk + i)^/ , , 1+iyk N 
c(k,q) = ^ -i-F ( + n,;1 - W — ^- | + 



(yk-i^ 



r(i-n) 

Now, according to [35J 



F -p., 1 + m., ; 1 - n; 



2 

1 — iyk 



F(-p, 1 + ml - n;z) = (1 - z)^F(1, -2 W 1 - mz), 



and then 



c(k,q) 



r(i - 



f i , —2m i - n, 



f( i,-2 W i-m ^-^| -i- 



1 — iyk 



so that via the result [36] we have 



F 1 , -2\i\ 1 - n; 



1 +iyk 



2F(1-^y^ 



FN, -2p.; 1 - |x; 



(1 +y 2 k 2 )^. 



1 — iyk 



Eqs. floTTTj) and fl5TT2|) lead to the cut f(k, q) of F(k, q) in (ET7J) 



f(k,q)= 1 



4cc 
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which is the q-Fourier transform of the q-Gaussian on the real axis for fixed 
q. This result can also be obtained using the real q-Fourier transform 



f(lc,q; 



C q [l + (q-1)cxx 2 ] — X 



1 + (1 - q)ikx{ C q [1 + (q - 1)<xx 2 ]^ 



q-1 1 i-q 



dx, (5.14) 



1 < q < 2. 



Conclusions 

By recourse to tempered ultra-distributions we have evaluated a complex q- 
Fourier transform F(k, q) for the q-Gaussian that is one-to-one, solving thus 
a flaw of the original F q — definition, i.e., not being of the essential one-to-one 
nature, as illustrated, for instance, in [25J. The computation has been done 
both for a floating q and for a fixed one. An essential piece of the Tsallis' 
machinery has thus been rebuilt. 
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6 Appendix: Tempered Ultradistributions and 
Distributions of Exponential Type 

For the benefit of the reader we give a brief summary of the main properties 
of distributions of exponential type and tempered ultra-distributions. 
Notations. The notations are almost textually taken from Ref. [21] ■ Let 
]R n (res. C n ) be the real (resp. complex) n-dimensional space whose points 
are denoted by x = (xi,x 2 , ...,x n ) (resp z = (zi,z 2 , ...,ZnJ). We shall use the 
notations: 

(a) x + y = (xi +yi,x 2 + y2, —,*n + TJn) ; ax = (axi, ax 2 , ax n ) 

(b) x ^ means Xi ^ 0,x 2 ^ 0, ...,x n ^ 

n 

(c) x -y=Y. XjUj 

j=i 

n 

(d) | x |= | xj | 

H 

Let N n be the set of n-tuples of natural numbers. If p G N n , then p = 
(pi , p 2 , p n ), and pj is a natural number, 1 5= j 5= n. p + q stands for 
(p! + qi,p 2 + q 2 , ...,p n +qj and p ^ q means pi ^ q b p 2 ^ q 2 , ...,p n ^ q n . 

n 

x p entails x v ^x\ 2 ...x^ . We shall denote by | p |= Y.V) an d ca ^ L) p the 
differential operator 3 p,+P2+ - +Pn /3x 1 Pl 3x 2 P2 ...3x n Pn 

For any natural k we define x k = x\x\...x\ and 3 k /3* k = d nk /dx\dx\...dx k 
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The space H of test functions such that e p|x| |D q (b(x)| is bounded for any p 
and q, being defined [see Ref. ([31])] by means of the countably set of norms 

||$|| = sup eP M |D<$(x)| , p = 0,1,2,... (6.1) 

0<q<p,x 

The space of continuous linear functionals defined on H is the space of 
the distributions of the exponential type given by ( ref. [21] )• 

T = w t ekWf (x) ] (6 - 2) 

where k is an integer such that k ^ and f (x) is a bounded continuous 
function. In addition we have 1-L C «S C <S' C A^, where «S is the Schwartz 
space of rapidly decreasing test functions (ref [32] ) . 
The Fourier transform of a function § G H is 

oo 

1 



2n 



4)(x) e tz ' x dx (6.3) 



According to ref. [31] , 4>(z) is entire analytic and rapidly decreasing on straight 
lines parallel to the real axis. We shall call H the set of all such functions. 

U = 7{U] (6.4) 

The topology in H is defined by the countable set of semi-norms: 



k = sup |z|>(z)|, (6.5) 

z6V k 
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where V k = {z = (z!,z 2 , ...,Zn) £ C n :| ImZj |^ k, 1 ^ j ^ n} 

The dual of H is the space U of tempered ultra-distributions [see Ref. ((31] )]. 

In other words, a tempered ultra-distribution is a continuous linear functional 

defined on the space "H of entire functions rapidly decreasing on straight lines 

parallel to the real axis. Moreover, we have M. C S C S CU. 

tl can also be characterized in the following way [see Ref. ((21] )]: let A. w 

be the space of all functions F(z) such that: 

A) - F(z) is analytic for {z G C n : |Im(zi)| > p, |Im(z2)| > p, |Im(Zn)| > p}. 

B) - F(z)/z p is bounded continuous in {z G C n : |Im(zi)| ^ p,|Im(z2)| ^ 
p, |lTa(Zn)| ^ p}, where p = 0, 1,2, ... depends on F(z). 

Let TT be the set of all z-dependent pseudo-polynomials, z G C n . Then W is 
the quotient space 

C) - U = A w /U 

By a pseudo-polynomial we understand a function of z of the form 
^ s Z-G(zi,...,Zj_ 1 ,Zj+i,...,Z T with G(zi,...,Z j _ 1 ,Zj+i,...,Z a ) G Ao 
Due to these properties it is possible to represent any ultra-distribution as 
[see Ref. ((3TJ )] 

F((j>) =< F(z), 4>(z) >= oF(z)4>(z) dz (6.6) 

r 

r = P] U Yz U ...r n , where the path F, runs parallel to the real axis from — oo 
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to oo for Im(zj) > C, C > P and back from oo to — oo for Im(Zj) < — C, 
— C < —p. (r surrounds all the singularities of F(z)). 

Eq. (16. 6p will be our fundamental representation for a tempered ultra- 
distribution. Use is also made of the "Dirac formula" for ultra-distributions 



[see Ref. ([30])] 

Ffzl 



(2m) 



f(t) 



(t! -Z 1 )(t 2 -Z 2 ) 



(t n - Z n J 



dt 



(6.7) 



where the "density" f (t) is such that 



oF(z)cp(z) dz 



f(t)4)(t) dt. 



(6.8) 



While F(z) is analytic on V, the density f(t) is in general singular, so that 
the r.h.s. of (16. 8p should be interpreted in the sense of distribution theory. 
Another important property of the analytic representation is the fact that 
on P, F(z) is bounded by a power of z [3T] 

|F(z)| < C|z| p , (6.9) 

where C and p depend on F. 

The representation (16.61) implies that the addition of a pseudo-polynomial 
P(z) to F(z) does not alter the ultra-distribution: 

c>{F(z) +P(z)}4>(z) dz = oF(z)4>(z) dz + oP(z)4)(z) dz 
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However, 



oP(z)(j)(z) dz = 0. 



As P(z)cb(z) is entire analytic in some of the variables Zj (and rapidly de- 
creasing), we obtain: 



(){F(z) + P(z)}cMz) dz = oF(z)cMz) dz. 



(6.10) 
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